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ABSTRACT 



This paper characterizes the physical and kinematic properties of external massive star-forming regions in a sample of (U)LIRGs. We 
use high angular resolution ACS images from the Hubble Space Telescope (HST) in F435W (~ B) and F814W (~ /) bands, as well as 
Ho line emission maps obtained with integral field spectroscopy We find 31 external Hof-emitting (i.e., young star-forming) complexes 
in 1 1 (U)LIRGs. These complexes have in general similar sizes (from few hundreds of pc to about 2 kpc), luminosities (Mf^^sw < 
-10.65 and L(HQf)>10^^Lo), and metallicities (12 + log (0/H) ~ 8.5-8.7) to extragalactic giant Hii regions and TDG candidates found 
in less luminous mergers and compact groups of galaxies. We assess the mass content and the likelihood of survival as TDGs of 
the 22 complexes with simple structures in the HST images based on their photometric, structural, and kinematic properties. The 
dynamical tracers used (radius-cr and luminosity-cr diagrams) indicate that most of the complexes might be self-gravitating entities. 
The resistance to forces from the parent galaxy is studied by considering the tidal mass of the candidate and its relative velocity with 
respect to the parent galaxy. After combining the results of previous studies of TDG searches in ULIRGs a total of 9 complexes satisfy 
most of the applied criteria and thus show a high-medium or high likelihood of survival, their total mass likely being compatible with 
that of dwarf galaxies. They are defined as TDG candidates. We propose that they probably formed more often during the early phases 
of the interaction. Combining all data for complexes with IFS data where a significant fraction of the system is covered, we infer a 
TDG production rate of 0.3 candidates with the highest probabilities of survival per system for the (U)LIRGs class. This rate, though, 
might decrease to 0.1 after the systems in (U)LIRGs have evolved for 10 Gyr, for long-lived TDGs, which would imply that no more 
than 5-10% of the overall dwarf population could be of tidal origin. 

Key words, galaxies interactions - galaxies: dwarf - galaxies: formation - galaxies: evolution - methods: observational 



1. Introduction 

The question of how galaxies form has been a key issue in 
extragalactic astronomy for the past few decades. One product 
of galaxy formation observable primarily in the local universe 
is the dwarf galaxy population, that usually forms in tidal tails 
during interactions of giant galaxies up to the present day. 
These galaxies are now generally called tidal dwarf galaxies 
(TDGs; Due & Mirabel 1998) . 

The scenario of low-mass galaxy format ion durin g gian t 
galaxy collisions was first proposed by Zwick^ (Il956h 
despite a lack of strong qbservat ional evidence, which was 

Subsequent observations 
1993) and numerical sim- 



^ . provided later (ISchweizeil Il978j) 



! (iMirabel et a L* "1991"; "Hernauid 



ulations dBarnes & Hernquist 119921: lElmegreen et alJ Il993h 
demonstrated that the existence of dwarf galaxies as a result 
of a major interaction is possible. Since then and up until 
now, several observing campaigns have been launched in 
intera cting systems, such as those in NGC 2 782 (lYos hida et all 
119941) . Arp 105 (Due & Mirab^ 1 19941: iDuc et al. 1995) 
NGC 7252 (Hibbardetalj|i994D, NG CJggTjbuc & Mirabei 
119981: Ikisdon et al. 2QQQ), Arp 245 d Duc et a l.1 |2QQQ|) small 



sample s ( Weilbacher et al] I200QL I2QQ3 : Knier man et aD l2003l: 
'Monreal-Ibero et al. 2007), and Arp 305 (HancocketaljL2009"). 
Jjibbard & Mihos (1995) also developed a successful dynamical 
N-body model of NGC 7252. Some other investigations have 
tried to find TDGs specifically in compact galaxy groups 



(e.g., 'Hunsberger et al."'l996'; 'Iglesias-Pa ramo & Vflche3l200ll: 
iTem porin et al. 2003; Nishiura et al. 2002^^ 

All spectroscopic observations of TDGs have shown that 
the luminosity-metallicity correlation found for normal dwarf 
galaxies does not hold for TDGs. While dwarf galaxies have 
lower oxygen abundances at lower luminosities, TDGs have 
an approximately constant metallicity of around one-third of 
the solar value as determined from gaseous emission lines (see 
e.g. jDuc et al ] [200Ql:lw^acher et al.l[20Q3h . 

The identification of the observed TDG candidates as real 
galaxies is not straightforward. They are embedded in a tidal 
tail and, most probably, do not have massive dark matter 
halos. Tidal forces produced by the parent galaxy disturb their 
gravitational field, strong star formation might blow away the 
recently accreted gas, and som e of the TDG s may even fall back 
into the central merger ( Hi bbard & MihosI 11995). An accepted 
definition that tries to ensure that only the objects that we refer to 
as tidal dwarf galaxies deserve the classification of a "galaxy" is: 
A tidal dwarf galaxy is a self-gravitating entity of dwarf-galaxy 
mass built from the tidal material expelled during interactions 
(iDuc et al.ll200a:lWeilbacher & Fritze-v. Alvenslebenll2 001). 

Taking into account evaporat ion and fragmen tation pro- 
cesses as well as tidal disruption, iDuc et aL 1 (l2QQ4l) suggested 
that a total mass of as high as 10^ M© may be necessary for 
it to become a long-lived TDG. This is the typical mass of the 
giant HI accumulations observed near the tip of several long tidal 
tails. Less massive condensations may evolve, if they survive. 
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into objects more similar to globular clusters. However, this 
mass criterion is not well-established, and objects with a total 
mass of a few 10^ -a few 10^ M© are normally considered as 
TDG candidates. 

Interacting systems (and remnants) constitute the optimal 
environment for finding TDGs. Although many studies devoted 
to searching and characterizing TDGs in nearby interacting 
systems have been carried out in the past few decades, only 
a few searches for TDG candidates have been performed 
to date in the most energetic interacting systems in the 
local universe, luminous (LIRGs) and ultraluminous (ULIRGs) 
infrared galaxies. 



infrared 



The LIRGs and ULIRGs are objects with 
luminosities of lO^^L© < L^oi ~ Ljr[^ - lOOOyum] < 
and IQ^^L© < L//?r8 - lOOOyumCl < lO^^L©, respectively 
(ISanders & Mirabeil ll996'). The main source of th is luminosity 
in mo s t of them is t hought to be starburst acti vity (iGenzel et al.l 
119981: iFarrah et al.1 [2003; Y uan et al.l l201Ql) . These galaxies 
are rich in gas and dust, and more than 50% of LIRGs 
and most ULIRGs sho w signs of being involved in a major 



interaction/merger (e . g Surace et al 



Farrahet aLl l200ll: [Bushouse ~ al.) l2002l: lEvans et al.1 



Veilleux et al.ll2006h 



1998L 120001: ICui etal 



2001 



20021: 



On the basis of previous studies of nearby interacting 
galaxies, TDGs are normally identified as extreme star-forming 
regions in tidal tails. Owing to the high star-formation rates 
measured in (U)LIRGs and the interacting nature of many of 
them, we expect to find a significant number of candidates in 
this environment. Previous s tudies of TDGs in UL IRGs include: 
one of the Superantennae (iMirabel et al.l 1 19911) , whose TDG 
candidate turned out to be a background object ('Mirabel et al.' 
Il99ll) : that of fo ur ULIRGs, wh ere only three TDG candidates 
were identified (iMihos & BothunI 1 1998 ): and another of nine 
ULIRGs where, in five of them, a total of twe lve TDG candidates 
were identified based on their Hof-emission ( Monreal-Ibero et aP 
[2QQ7h . This third work also studied the likelihood that the 
TDG candidates will survive, and found that five out of the 
twelve show high probabilities of t heir remaining a TDGr. This 
paper extends the previous study of iMonreal-Ibero etaL I (l2007h 
by including galaxies in the LIRG class. 

A study of TDGs in the local universe is issential to aid our 
understanding of galaxy formation at high-z. If TDGs, which 
form from tidal debris, are long-lived, they could contribute 
significantly to the total population of dwarf satellites, in 
addition to primordial dwarfs. Their statistical properties would 
then have to be modified and the available constraints on 
cosmological models would have to be updated. Early studies 
showed that this is the case, e ven that the overall dwa rf 
population could be of tidal origin (lOkazaki & Taniguchi 2000). 
However, it has been claimed that only a marginal fraction 
(less than 10%) of dwarf galaxies in the local universe could 
have a tidal origin (iBournaud & Dud l2006l: IWen et al.l l201ll: 
lKaviraietal.ll201lh . The (U)LIRGs are major co ntributors to 
the st ar-formation rate density at z~l-2 (Perez-Gonzalez et al. 
l2005h . In addition, the observed properties of (U)LIRGs share 
many similaritie s with populatio ns of star-forming g alaxies at 
higher redsh ifts (|QiamnaiLeLayi2QQ3|; iFraver et al.ll2003L .2004: 
lEngel et al.l l2010|). Hence, the study of the local (U)LIRG 
population provides an opportunity to link the properties of 
high-z galaxies with those we observe in the nearby universe. 



^ For simplicity, we define the infrared luminosity to be 
LiR (Lo ) = log (L;^[8 - mOfim]) . 



In particular, they can be the key to understanding how dwarf 
galaxies form in the early Universe. 

In this paper, we characterize the extranuclear star-forming 
regions of a sample of local (z < 0.1) (U)LIRGs using 
information derived from integral field spectroscopy (IFS) data 
together with that of high resolution images from the Hubble 
Space Telescope (HST) in the B and / photometric band. We 
perform a dedicated search for potential TDG candidates in local 
(U)LIRGs. We derive and compare candidate properties such as 
the metallicity, age, and mass of the young stellar population. 
Ha luminosity, dynamical mass, velocity dispersion, relative 
velocity etc., with those of TDG candidates in the literature. We 
use these properties to estimate the likelihood of the survival of 
these regions as future TDGs. 

The paper is organized as follows. We first describe the 
sample in section [2l The photometric and spectroscopic data 
are presented in section [S] We then identify and characterize 
the luminosities, colors, metallicities, etc. of extranuclear Ha- 
emitting complexes in IH In section [3 we select the complexes 
that are most likely to constitute TDG candidates and discuss 
their total mass content and the likeliness of the survival of these 
regions as future TDGs, and the implications for the formation of 
TDGs at high-z. Finally, we draw our conclusions in section [6l 

2. The sample 

The galaxies selected for this study were taken from a 
representative sample of 32 low-z both luminous and 
ultr aluminous infrared g alaxies, (U)LIRGs, presented 
in iMiralles-Caballero et al.l (12011) (hereafter MCll). It 
covers the luminosity range 11.39 < Ljr < 12.54 and the 
distance range from 65 to about 550 Mpc. It also spans all types 
of nuclear activity, with diff'erent excitation mechanisms such 
as LINER (i.e., shocks, strong winds, weak AGN), HII (star 
formation), and Seyfert-like activity (presence of an AGN). 
Finally, all the diff'erent morphologies usually identified in these 
systems are also sampled. 

The main purpose of the paper is to combine photometric 
and spectroscopic data to perform a physical characterization 
of the most luminous extranuclear Hof-emitting regions in these 
galaxies. Therefore, we selected systems for which IFS data was 
available at the moment of analysis, that is a total of 27 systems. 
As explained in section 14. 1[ regions of interest were selected 
as any high surface brightness compact region (Ha clumps) 
in the emission line maps (obtained from the IFS data) at a 
projected distance from the nucleus of the galaxy greater than 
2.5 kpc, and associated with one or several of the star-forming 
regions (knots) in the images taken with the Advanced Camera 
for Surveys (ACS), onboard the Hubble Space Telescope (HST). 
On the basis of these criteria Hof-emitting complexes have been 
identified in 1 1 systems. We did not identify other complexes 
in the remaining sample basically for two reasons: (1) based on 
the above criteria, in some systems knots are detected but no 
Ha emission is found in the outskirts in the Ha line map (e.g., in 
IRAS 22491-1808); (2) the field of view (FoV) covered by the 
IFS data is not large enough to encompass the full extent of some 
systems. Hence, the nature of blue knots along the tidal tail and 
at its tip, where TDGs are normally found, cannot be assessed 
using the IFS data. 

Therefore, the final sample consists of 11 systems with at 
least one region of interest (see table [TJ, or more precisely 7 
LIRGs and 4 ULIRGs. For completeness, we also include in 
our tables the results of iMonreal^ Ibero et al.l (l2007h for IRAS 
16007-1-3743, for which an image in the B photometric band was 
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unavailable. According to the classification scheme defined in 
MCll, IRAS 04315-0840, IRAS 08355-4944, IRAS F10038- 
3338, and IRAS 15250-h3609 are in an advanced phase of 
a merger (merger phase) where the nuclei have apparently 
coalesced. However, the tails are still recognizable in the images, 
in some cases very prominently. The rest of the systems are in 
earlier stages of the merging process, first contact of pre-merger 
phases, according to the classification scheme we use. 

3. The data 

The results presented in this paper are based on photometric and 
spectroscopic data. In the following, we give the details of both 
sets. 

3.1. HST imaging 

We retrieved high angular resolution archival HST images taken 
with ACS in two broad-band filters, F814W and F435W. The 
former is equivalent to the ground-based Johnson-Cousins / 
and the latter diff'ers from the g round-based Johnson -Cousins B 
between 7% and 12% in flux dSirianni etal]l2005h . The pixel 
size for this instrument is 0.05'', and given the FoV of the 
instrument (200'' x 200") the galaxies and their tails in the / 
photometric band are completely covered. More details about 
the observations and the post-calibration applied can be found in 
MCll. 

Complementary images from the Near Infrared Camera and 
Multi Object Spectrometer (NICMOS), onboard if^T, were also 
retrieved. They were taken with the F160W (H) band filter. For 
IRAS 06076-2139, we used an available H band image taken 
with the Wide Field Camera 3 (WFC3). These H images helped 
us determine eff'ective radii of the galaxies. 

3.2. Integral field spectroscopy maps 

We obtained IFS observations for the galaxies of the sample with 
two diff'erent fiber-based optical integral field systems: 

- The ULIRGs in our sample were observed with 
INTEGRAL ( Arribas et'aP [19 98). connected to 
WYFFOS (Bingham eL^niaal and mounted on the 
4.2 m William Herschel Telescope. The FoV varies 
depending on the INTEGRAL dithering and pointings used, 
though in general it is around 16"x 12.3". The angular 
sampling is 0.90"per fiber. Details of the data r e duction and 
calibration can be found in iGarcia-Marin et al.l (l2009bh and 
references therein. 

- The LI RGs in our sample were observed with 
VIMOS (iLeFevre et al.l |^03), mounted on the Nasmyth 
focus B on the VLT. The FoV covers 27" x 27" with an 
angular sampling of 0.67 "per fiber. Details of the data 
reduct ion an d calibration can be found in both lArribas et aP 
(120081) and lRodriguez-Zaurin et all (l201lh . 

In both sets of data, the Ha line lies within the spectral 
range and, once calibrated. Ha line maps were generated 
by fitting the line using Gaussian functions. Maps with the 
Ha equivalent width and the [Nii] lines were also generated. 
In addition, maps with the oxygen ([Oiii] /1/14959,5007), and 
HJ3 lines for galaxies observed with INTEGRAL were produced. 
In this paper , we use the maps publi s hed in Garcia-Marm et al. 
([2009b) and iRodriguez-Zaurin et al.l (l201lk for the systems in 
our sample. 



4. Results 

4.1. Identification of the Ha -emitting complexes 

The identification of the stellar population that ionizes the 
surrounding interstellar medium, and subsequently produces the 
Ha emission, is not necessarily straightforward. Depending on 
the resolution of the data and the contamination of the field by 
other ionizing sources, the establishment of this matching is not 
always possible. We first selected bright condensations in the 
Ha maps that were coincident with at least one blue knot in 
the F435W ACS images. We chose to use these high resolution 
images instead of the F814W ones because the most likely 
responsible for the ionization were detected most efficiently with 
the blue filter. 

We defined a complex of young star formation as an 
Ha clump (a set of spaxel^ that had coincident position with 
one or several knots identified in the HST images. In general, 
around 12 spaxels (e.g., 3 x 4 spaxels) are needed to define a 
complex, depending on how isolated it is and the extent of the 
Ha emission before it is aff'ected by another source or diff'use 
nuclear emission. In one case (IRAS 08572N) a prominent 
Ha peak is observed at around 7 kpc west of the northern 
nucleus, but with no corresponding region in the HST images. 

We visually inspected each Ha map and found 1 1 systems 
with prominent peaks apart from the nuclei (see Fig. [T]). We 
defined the inner complexes as those located within a projected 
distance of 2.5 kpc to the closest nucleus and the outer 
complexes as those further out than 2.5 kpc. Given the resolution 
of the spectral maps (a factor of more than a hundred in area), 
the identification of the ionizing stellar knots is highly uncertain 
in inner complexes, since the Ha flux there is very much likely 
to be contaminated by flux from nearby knots. Together with the 
motivation of searching for TDG candidates, they constitute the 
main reason why we exclude these complexes from the study. 
Therefore, we focus on the outer complexes. In addition, we 
do not select complexes in rings of star formation (i.e., those 
at the external ring in IRAS 06076-2139), since they are not 
normally associated with tidal forces, hence with TDGs. All 
considered, we identified for our study 31 outer young. Ha - 
emitting complexes of star formation, shown in Fig.[T]and[2l 

4.2. Structure and location of the complexes 

Owing to the higher angular resolution of the HST images, we 
identified the Ha clumps (observed in the spectral maps) with 
a complex of star-forming knots (observed in the photometric 
images). The number of knots in the ACS image per Ha complex 
is on average 2.2. However, the structure of these clumps is in 
general quite simple. As can be seen in Fig. [2] and table [TJ the 
position of the Ha-emission peak is coincident or quite close 
to a prominent blue knot. In most cases, only one bright knot 
is observed within the area defined for the Ha complex (e.g., 
CI in IRAS 08572-h3915 N) or a centered bright knot plus a 
few fainter knots located at the border of the Ha complex (e.g., 
CI in IRAS 04315-0840). In a few cases, the structure is quite 
complex, with several bright blue knots spread inside the area of 
the clump (e.g., C2 in IRAS 12112-h0305). 

When characterizing the Ha complexes, given the rich 
structure of some of them, we assume that every knot detected in 
the HST images within the area of the complex is an Ha emitter. 
With the current data, the estimate of the contribution of each 

^ A spaxel corresponds to the spatial element from which we have 
one independent spectrum, and it is used to reconstruct the map. 
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Fig. 1: Systems under study with the Ha maps obtained from integral field spectroscopy. The F435W ACS images are shown and contours corresponding to the Ha are over- 
plotted. The contours are in arbitrary units, optimized to show and identify the structures of the Ha complexes. In some cases, there are two diff'erent pointings for the same system. 
The labels show the location of the complexes under study. The horizontal line at the bottom left corner shows a scale of 2.5 kpc. The VIMOS FoV is shown approximately. 
For INTEGRAL sources, the INTEGRAL FoV is over-plotted in the blue boxes. Most of the Ha peaks match up with several knots identified in the blue images. Note that the 
brightest knot on the Ha clump below CI corresponds to a red foreground star, not a knot of star formation. 



IRAS 04315-0840 IRAS 06076-2139 IRAS 07027-6011 S IRAS 08355-4944 IRAS 08572+3915 N 
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Fig. 2: Zoomed view of all the identified Ha complexes. In each case, the FoV has been selected to visualize the complex more clearly. The horizontal line at the bottom left 
corner shows a scale of 2 kpc (blue) for complexes detected on INTEGRAL maps and 1 kpc (white) for complexes detected on VIMOS maps. 
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knot to the total Ha emission or if a knot is not an Ha emitter 
(i.e., it is older than 10 Myr, at which the Ha emission declines 
abruptly) for all the complexes, is not possible. We also assume 
that only the detected knots are responsible for the whole 
Ha emission of the complex. Throughout the paper, we analize 
the validity of this hypothesis, since a priori knots below our 
detection limit in the ACS images might also contribute to the 
Ha emission. 

The Ha complexes are located at distances in the range 
3.3-13.7 kpc from the nucleus of the parent galaxy, with a 
median value of 9.3 kpc (see table [T]). These distances, owing 
to projection effects, represent lower limits to the real distances. 
They are relatively close distances to the parent g alaxy compared 
to those of the TDG cand i dates(35-100kpc; e.g ., lDuc & Mirabell 
ll998l:ISheen et al.l l2009l: iHancock et all l2009l) . although some 
TDG candidates have been reported at such distance s 
(e.g.. Iglesi as-Paramo & Vflchezll2QQl]:IWeilbacher et al.ll2003l) . 
At these distances, some of the Ha complexes lie along the 
tidal tails (see Fig. [T]), and the diffuse low surface brightness 
Ha emission (which is not subtracted in our measurements of 
the Ha flux) is either undetected or considerably diminished, 
thus does not contaminate the emission from the knots. Others 
also appear to be simply a non-nuclear giant Hii region. 



4.3. Characterization of the Ha -emitting complexes 

We now investigate the physical properties of the bright 
Ha complexes identified. This characterization includes the 
analysis of the stellar continuum magnitudes (Mf435w and 
color), spectral features (Ha luminosity and equivalent width), 
as well as physical sizes and an estimate of the metallicity. 
All these observables are compared to those measured in 
extragalactic Hii regions, dwarf galaxies, and other TDG 
candidates. 



4.3.1 . Photometric properties of the complexes 

4.3.1.1. Broad-band luminosities and colors 

The integrated blue absolute magnitude Mf435w and color 
(Mf435w-Mfsuw ) of the young stellar population within 
the Ha complexes were estimated by adding all the 
flux of the knots inside a complex. We cover a mag- 
nitude range Mf435w = [-9.32,-15.77], and a color range 
Mf435w-Mfu4w = [-0.23,2.34] (see tabled]). The median values 
correspond to -12.06 and 0.63 magnitudes, respectively. The 
typical uncertainties are between 0.05 and 0.1 magnitudes. 

The integrated luminosities are in general higher than in 
embedded clusters in nearby extragalactic Hii regions in spirals 
(iBresolin & Kennicutill997l) and more typical of those found 
in ex tragalactic giant Hii regions ('Ma yyalll994l: iFerreiro et al.l 
l2008l) . Although early studies in extragalactic Hii regions only 
focus on nearby galaxies (Z)l < 20 Mpc), Ferreiro and co- 
workers observed minor mergers at distances (Z)l = 40-170 
Mpc) similar to our sample, detecting also only the most 
luminous giant Hii regions, the external ones covering the range 
Mb ^[-11.7,-17.4], typically 2 magnitudes more luminous 
than the range sampled in this study. 

Nearby dwarf galaxies are typically more 
luminous than B < -13.0, with colors in th e range 
MF435w-MFf^u w = 0.7- 2 mag ([H unter & Gallaghe j 1 19861: 
iMarlowe et aP 1 19971: ICairos et a l. 2001 ). In a d dition 
the luminositi e s ob se rved for nearb y (iDuc et al.l 119971: 
iDuc & Mirabell 119981: iDuc et al.l l2007h and more distant 



(IWeilbacher et al.ll2000l: pTemporin et al.ll2003l) TDG candidates 
are generally more luminous than B = -10.65, as for most of our 
Ha complexes. 

The integrated broad-band luminosities of the embedded 
stellar population in our Ha complexes is, therefore, compatible 
with the luminosities measured in giant Hii regions and TDG 
candidates. 

4.3.1.2. Sizes and compactnesses 

Owing to the higher angular resolution of the HST images, 
estimates of the sizes and compactnesses of the Ha emitting 
clumps were derived from the ACS images. All the knots 
within a complex were assumed to represent the ionizing 
stellar population. We then added up the area of each knot 
(derived in MCll) within an Ha complex, computing a total 
area At. With this a rea, we derived an equivalent total radius, 
given by r= ^JA^|7l (see table [T]). Had we used the spectral 
maps to measure the radius of the complexes {rwa in table [T]), 
we would have overestimated it in general by more than a 
factor of two. The radii r range from somewhat less than 100 
pc to about 900 pc, with a median value of 280 pc. These 
sizes are similar to those for the largest giant extra galactic 
Hii regior 
119941: iRozas et al 
ground-based instrumentation. 

To evaluate the compactness of the Ha complexes, we also 
derived an effective radius (reff) for each, as an approximation 
of the half-light radius. We again used the HST images for 
this determination. In this case, since in most cases one knot 
dominates the total broad-band luminosity within the complex, 
we identified the effective radius of an Ha complex with the 
effective radius of the most luminous knot (derived in MCI 1). In 
complexes where more than one knot dominates the luminosity, 
this approach lead to an underestimate of the effective radius. 

The effective radii of complexes observed with VIMOS 
are smaller than 100 pc (see table [TJ, from 20 pc to about 
80 pc. For those observed with INTEGRAL, the values range 
from 100 pc to about 300 pc. These sizes are in general larger 
than the effec tive radii of the so-c alled ultra compact dwarf 



Hii regions observed in nearby (Z)l < 20 Mpc) spirals (iMayyal 
, '2006a'; 'Dia z et al.l |2007|), measured using 



galaxies (e.g.. iDabringhausen et al.h2008) , similar to those in 
giant Hii regions fe.g lRelafio et al] 120051: iRozas et al.l l2006ah 
and within the range of the smallest blue compact dwarf galaxies 
(with Teff =0.2-1.8 kpc;e.g., Cairos et al. 2003; Papaderos et alj 
2006; Amor m et al.l 120091) and Local, N GC 1407, and L eo 
groups (with reff ~ 0.3-1 kpc: lMateollT998HFbrbes et al.1l201lh . 

The ratio of the effective to equivalent radii gives an idea of 
the compactness of a certain complex. They are very compact, 
with a ratio of typically reff/r = 0.2, ranging from 0.15 to 0.6. 



4.3.2. Ha luminosities and equivalent widths 

Both Ha flux and equivalent width measurements (EW) 
were obtained for the 31 Ha complexes. Their observed 
Ha luminosity is typically higher than 10-^^ erg s"^ (see 
table 12]) . The luminosity of about 20% of them is even 
higher than 10"^^ erg s"^ The median value of the Ha flux is 
2.9x10^^ erg s"^ Monr eal-Ibero et al.l (l2007b also measured the 
Ha fluxes for four of our systems with INTEGRAL data (in 
IRAS 08572-h3915 only for the northern pointing) and most of 
them are systematically lower by up to a half. The discrepancy 
occurs because they used apertures of 0.5''or 1.0'', whereas 
we have added up the fluxes of several spaxels, considering 
a greater effective aperture (one spaxel on INTEGRAL maps 
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Fig. 3: Relation between the Ha luminosity and the size 
of the identified Ha complexes in our sample (black dots), 
extremel y luminous Ha co mplexes in the Antennae (blue 
plusses; Bastia n et alJ |2006), diff'erent samples of isolated 
giant Hii regions in spirals and minor mergers (red sym- 
bols; iTerleyich & MelnicS 1 19811: iFerreiro & Pastorizal I2QQ4I: 
iRozas et alJl 2006a'). TDG candidates in CO J1720-67.8 (green 
diamonds; [lemporin et al. 2003) and nearby dwarf amorphous 
(yellow crosses; Marlowe et al. 1997). None on the luminosities 
on the plot are corrected for internal extinction. 



already corresponds to an aperture of 0.45'' in radius). In the 
case of IRAS 121 12+0305, we grouped their kc and kl regions 
and did not consider R2 because it is too close to the northern 
nucleus, which makes its Ha flux determination very uncertain. 

Fig. [3] compares the Ha luminosities and sizes of our 
complexes with those of extragalactic Hii regions. Ha complexes 
in nearby galaxies, TDG candidates, and dwarf galaxies. This 
plot shows that using the Ha luminosity alone to establish the 
nature of a given object might be misleading. For instance, 
the complexes of star formation in the Antennae have an 
Ha luminosity comparable to the brightest extragalactic Hii or 
even to dwarf galaxies, while these complexes are generally 
much smaller. Our complexes are typically located where the 
giant extragalactic Hii regions and the dwarf-type objects lie on 
the plot. 

We have also estimated the total equivalent width of a 
given complex from the EW and the Ha spectroscopic line 
maps. Since in a spaxel the EW oc Flux Ha / Flux continuum, 
we first determined the continuum per spaxel. We then added 
up all the continuum within a complex to compute the total 
continuum of the complex. Finally, we divided the total Ha flux 
by the total continuum flux for each complex. Their integrated 
EW span a range of two orders of magnitudes, from 4.1 A 
to about 170 A (see table O. For each complex, we also give 
the value of the spaxel with the highest EW (EWpeak). The 
highest value corresponds to EWpeak=374A. Its median value 
(< EWpeak > = 85A) doubles the median value of the integrated 
EW computed for the complexes. The Ha equivalent widths 
derived, which is indicative of a very young stellar population, 
are comparable to those measured for extragalactic Hii regions 
and TDG candidates (see Fig. [4]). 

The ratio ew\h^ defined for each complex, where 
PeakEW refers to EWpeak- This ratio helps us understand how 



Fig. 4: Equivalent widths measured for the com plexes in this 



study (black), for extragalacti c Hn regions (green :lMavvd 1994 ) 



and for TDG candi dates (red; llglesias-Paramo & Vflchezll2001 



iTemporin et al.ll2003D .The number of objects is indicated in each 



case. 



the ionizing population is distributed. If two complexes have 
a similar underlying old population, the ionizing population in 
the one with the larger ratio is probably more concentrated (i.e., 
less knots) than the other with the smaller ratio. And if it is 
less concentrated (i.e., has more knots), most of the knots are 
probably somewhat older (low EW) and only a few are very 
young (high EW), leading to an older average population. In fact, 
the complexes with large (> 3) ratios have typically few knots 
(1-3), such as CI in IRAS 06076-2139, CI in IRAS 08572+3915 
N, and C4 in IRAS 12112+0305. 



4.3.3. Contamination of the Ha luminosity from embedded 
young populations 

With the knowledge of the integrated broad-band luminosities, 
the Ha luminosities, and the equivalent widths of the complexes, 
we can try to evaluate whether the measured Ha flux is emitted 
by possibly embedded (and undetected) star-forming knots. The 
existence of undetected Ha -emitters is justified with the Ha- 
emitter in IRAS 08572+3915 N (clump in the north-west). 

We estimated how many undetected knots there can be with 
an age derived from the peak values of the EW (typically 
between 5 and 10 Myr, according to the stellar population 
models used in section [SXT]) that can contribute significantly to 
the total broad-band luminosity of the knots. Since they would 
have a similar age, an embedded population of knots with a 
combined broad-band flux similar to that of the detected knots 
would emit half of the Ha flux. We then estimated how many 
knots at the detection limit would be needed to double either 
the F435W or the F814W flux. A total of 9 and 75 undetected 
knots per complex would be needed and we detected an average 
of 2.2 knots per complex. Thus, between 4 and 37 times more 
knots fainter by 2.4-4.7 magnitudes would be needed only in 
one complex. Assuming a luminosity function with slope of 2 
we would need practically one-third of the total knots predicted 
for the whole galaxy inside a complex, which is unlikely. 
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Fig. 5: Metallicity-luminosity relation for external 
Hii regio ns in Stephan' s Quintet (yellow open 
circles; Mendes de Oliveira et al. 2004), nearby irregular 
galaxies (b lack plusses, Pilvusin et al. 2004; black 
asterisks, van Zee & HaynesI 120061). TDG can didates 
(gree n filled tri angles. IWeilbacher et aP |2000 !. 
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squares, Tempo rin et al 



Due & Mirabell 1 19981: blue 
2003h and our complexes 



dots with errors). The line shows the correlation found by van 
Zee. Only the luminosities given by Pilyugin are corrected 
from internal extinction. For those Ha complexes for which an 
estimate of the extinction has been made (see section [SXTI) . the 
metallicity-luminosity relation is also shown (dots with dashed 
line errors). 



4.3.4. Metallicities 

Numerical studies have claimed that during an interaction large 
amounts of gas flow toward the central regions, carrying less 
enriched gas from th e outs kirts of the galaxy int o the central 
regions (iRupke et al.l 1201 Ol: iMontuori et al.l l2010l) . Indications 
of this mixing process, which usually disrupts any metallicity 
gradient and dilutes the ce ntral metall i city, h ave been observed 
in merging (U)LIRGs (Rupke et al. 2008 ), in galaxy pairs 
dKew lev et al. 2006; Ellison etal. 2008; Michel-Dan sac et al.l 
l2008h . and in star-forming galaxies at high-z (ICresci et 'al. 2010^ 
As a consequence of the metal enrichment of the external 
regions of the mergers, additional star formation can occur 
from reprocessed material there. We investigated whether the 
interaction process afl'ects the metallicity of our Ha complexes. 

Abundances are usually estimated using empirical methods 
based on the intensities of several op tical lines. The most 
popul ar methods are the widely used R23 (iTorres-Peimbert et"aD 
119891) and S23 ( Vflchez & Estebanlll99 6'') calibrators. However, 
both methods use emission lines that our spectral range 
does not cover. We us e d ins tead the N2 calibrator pro- 
posed by 'Denicolo et al.' J2002h and the empirical diagrams 
offEdmunds & Pagel (198%. The former is based on the ratio of 
the [Nii]/16584 to the Ha emission lines and the latter relates the 
ratio [Oiii]/14959 4 - rOniUSOO? toHj3 to the oxygen abundance, 
parametrized as in lDuc & Mirabell (Il998h . 

Given the spectral range of the VIMOS data, metallicity 
determinations have only been possible using the N2 calibrator 
for systems observed with this instrument. The typical 
uncertainties in the determined metallicities are about 0.2 dex. 



For the candidates observed with INTEGRAL, given these 
uncertainties, the difl'erences between both indicators (usually 
within 0.2 dex) are irrelevant, thus we adopted the average 
value. The adopted metallicities for all the candidates are 
shown in the last column in table [2l With values compatible 
to a solar metallicity, they generally range from 12 -h logO/H 
= 8.5 to 8.8 (with the exception of two candidates with 
lower abundances). The values for complexes observed with 
VIMOS have to be considered with care and must be checked 
with other metallicity indicators, since the N2 calibrator 
might be affected by ionization fr om strong shocks in th e 
external regions of these galaxies (Mon real-Ibero et al]|2010l) . 
In addition, spatially resolved studies of star-forming regions 
indicate that the assumption of spherical geometry is unrealistic 
in most cases, which has a d irect impact on the derivation 
of metallicities. lErcolano et al.l (2007) estimated the systematic 
errors in the metallicity determinations when assuming spherical 
geometry using difl'erent calibrations. In the worst-case scenario, 
the derived oxygen abundances might be overestimated by 0.2- 
0.3 dex by the use of the calibrators considered in this paper, in 
which case the derived metallicites for the Ha complexes would 
then approach Zo/3. 

The Ha complexes in this study do not follow the well- 
known metallicity-luminosity relation for nearby isolated dwarf 
galaxies (see Fig. O. Most galaxies also follow this relation 
( Weilbacher et al. 2003). On the other hand, Hii regions in 
compact groups of galaxies and TDGs in general deviate 
significantly from the relation and have a metallicity that is 
independent of their luminosity, an indication that all these 
objects consist of recycled material. 

5. Discussion 

Most Ha -emitting complexes have similar observational 
properties (i.e.. Ha equivalent widths and luminosities, 
Mf435w magnitudes, metallicities, radii) to the most luminous 
extragalactic giant Hii regions in spirals and mergers, as well 
as more massive objects such as TDGs or TDG progenitors. 
Associations of young star-forming regions with a larg e HI 
reservoir have been found in TDGs, for instance by Due et al.l 
( 2007). We may then consider whether the Ha -emitting 
complexes in our sample are dynamically unbound associations 
of objects with masses similar to observed super star clusters 
(SSCs) or larger, more massive, and self-gravitating objects such 
as dwarf galaxies. 

TDGs are self-gravitating objects with masses and sizes 
typical of dwarf galaxies (i.e., a total mass of 10^-10^ M©), 
formed with recycled material from the parent galaxies involved 
in an interaction/merger (Due et al. 2000). Hence, they are 
stable entities with their own established dynamical structure. 
To evaluate whether our complexes constitute real TDGs or 
TDG progenitors and assess their chance of survival, we need 
to answer a few basic questions: Is the complex massive enough 
to be considered as a dwarf galaxy? Is it stable enough to be 
unaff'ected by its internal motion? Are the gravitational forces 
too strong to disrupt it? In this section, we disscuss these 
questions. 



5.1. Selection of Ha-emitting likely to represent TDG 
candidates 



Although old TDGs have been observed (^e.g.. lDuc et al ] l2007h . 
their identification is not straightforward. When an old TDG is 
detected, the tail from its place of origin will probably have 
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completely disappeared and the TDG may be classified as 
another type of dwarf galaxy. A determination of the metallicity 
and, especially, the total mass of the objects is needed to 
establish its tidal origin and mass. With spectroscopic data, we 
are able to study (though with some limitations and biases) the 
metal content, the dynamical mass, and diff'erent methods to 
establish whether a given candidate can withstand both internal 
and external forces. It is quite complicated to identify luminous 
enough condensations of old populations (i.e., very massive) and 
obtain kinematic information because the stellar absorption lines 
that would have to be used often have too low a signal-to-noise 
ratio (S/N). Although TDGs containing only old populations 
may be present in our sample, we are unable to detect them. 
Another way of searching for TDG candidates is based on 
either analyzing the Ha emission clumps (for that we have 
a sufficiently high S/N) or/and combining spectral and photo- 
metri c data of blue Hor-emitting objec t s in interacting systems 



(e.g.. Ilglesias-Paramo & VilchezI 120011: IWeilbacher et 



g systems 
^ 120031: 



iTemporin et al.l 120031: iLopez-Sanchez et al .1 12004 . Therefore, 
the selection of TDG candidates in our sample of (U)LIRGs 
is based on the detection of Ha emission clumps with 
stellar counterparts, which is similar to our identification of 
Ha complexes. 

However, we are limited by the angular resolution of the 
spectral data. To make our estimates as accurate as possible, we 
do not consider in our discussion here complexes with multiple 
luminous star-forming knots. Different types of motions within 
the components of a given complex can strongly affect the 
derivations mentioned before. The most characteristic examples 
of complicated complexes are CI in IRAS F10038-3338 and C2 
in IRAS 12112-h0305. This does not mean that the possibility 
of these complexes being TDGs is ruled out, just that we do not 
have enough resolution to study the kinematics of each knot, thus 
determine whether they are kinematically bound or not. 

We consider only complexes with simple structures, that 
is, mainly where only one knot is detected. We also include 
complexes with several knots, one of which (normally centered 
on the Ha peak emission) dominates the broad-band luminosity. 
The complex CI in IRAS 043 15-0840 would be a good example. 
On the basis of these criteria, the kinematic and dynamical 
properties of 22 complexes are derived and compared with those 
expected for a TDG. 

If we detect more than one knot within the selected 
complexes, we consider only the broad-band luminosity of the 
brightest knot. In these cases, the Ha luminosity of the complex 
is scaled in the same way as the F435W flux is scaled to the total 
F435W broad-band flux of the complex in table [T] In practice, 
we multiply the Ha luminosity by a factor that in general is 
higher than 0.5 (a factor of 1 means that there is only one knot 
emitting the whole broad-band flux, thus responsible for the 
whole Ha emission). The kinematic measurements, such as the 
velocity dispersion, are not performed for the whole complex, 
but only for the spaxels close to and located at the peak of the 
Ha emission. This permits us to avoid any contamination by a 
faint knot in the border of the complex (e.g., faint knots within 
CI in IRAS 04315-0840). 

Using the kinematic information provided by the IFU data 
as well as the photometric measurements, we assess the origin 
and likelihood of survival of the 22 selected Ha complexes, 
and evaluate whether they can be the progenitors of tidal dwarf 
galaxies. 



5.2. Are the selected Ha complexes massive enough to 
constitute a TDG? 

5.2.1 . Age and mass estimates of the young population 

The Ha emission can be used to constrain the properties of 
any recent episodes of star formation. Given the youth of 
the population we study, we use the si ngle stellar population 
synthesis model Starburs t99 (SB99, iLeitherer et al.1 1 19991: 
I Vazquez & Leithere3l2005h . which have been optimized for the 
analysis of young populations, to estimate its age and mass. We 
consider instantaneous burst models with a Kroupa IMF ( Kroup^ 
2002) over the range 0.1-120 M© and solar metallicity. These 
models are normalized to 10^ M©. Mass estimates are multiplied 
by a factor of 1.56 when the Salpeter IMF is considered. 

According to the SB99 models, the detection of an Ha flux 
and EW larger than about 13 A (see table [2] for specific values) 
implies that the stellar population is younger than 10 Myr. With 
the broad-band luminosities of the brightest knot within the 
complex and the scaled Ha luminosity, we can estimate to first 
order the age and the mass of the young population of the 
selected complexes. We note that, owing to the distance of the 
systems (Dl < 270 Mpc) and the shape of the HST filters, the 
Ha emission-line does not contaminate these broad-band filters. 

We computed the mass-independent evolutionary track 
shown in Fig. [6] using information about the ratio of the flux 
of the modeled population in the broad-band filters to the flux of 
the ionized gas in the Ha line. Placing these ratios on a plot for 
all observed complexes, we now estimate the age of the young 
population. Once the age is known (hence a M/L ratio), the mass 
can be derived. However, before going any further and based on 
the assumptions made, we first assess the sources of uncertainty 
that a priori could affect our results: 

- The internal extinction is unknown. In general, the 
extinction in the inner regions is patc hy with high peaks in 
(U)LI RGs , up to Ay ^ 8 mags (see lAlonso-Herrero et aP 
120061 and iGarcia-Marin et al.1 l2009ab . though it decreases 
considerably with the galactocentric distance. Were 
complexes to have low extinction values (Ay ^ 1.5mag) 
their flux quotients would tend to follow the evolutionary 
track (solid line in Fig. O, and the age would remain 
practically unchanged. Only one complex probably has a 
higher internal extinction (CI in IRAS 14349-1447). A high 
extinction value of Ay -3.5 mags was inde ed measured 
for this complex in Mo nreal-Ibero et al.l (120071) . The typical 
extinctions of the TDG candidates studied in Monreal-Ibero 
are between Ay = 0.7 and 2 mag, which is consistent with 
the extinctions derived for our selected complexes (see 
Fig.O. 

- The local background flux (from the underlying parent 
galaxy) assumed to be associated with either older 
populations or non-Hof line-emitters, was subtracted when 
the photometry was performed for the knots. However, a 
knot itself can be formed by a composite of young and 
old population. If we were to assume that a significant 
fraction of the red and blue fluxes measured for the knots 
originated from an older population, the evolutionary track 
would change (see the composite populations in Fig.O; the 
young population would then be even younger than initially 
predicted and consequently less massive. Almost half of the 
complexes are incompatible with the composite track, since 
their flux ratios would be too small, even smaller than for a 
zero-age young population. The other complexes, for which 
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Fig. 6: Flux ratios of the blue and red filters when the 
Ha emission can be measured in a complex. The broad-band 
filter fluxes have been multiplied by the band-width to obtain 
the values with the same units as the Ha flux. Solid line shows 
the evolutionary track of the same ratios for a single burst 
using SB 99 models (see text). Dotted line shows the track for 
a composite population where 99% of the mass corresponds to 
a IGyr population. The small triangle almost at the beginning 
(zero-age young population) of the second track indicates that 
the young population is 2Myr old and we add 1 Myr to the 
subsequent triangles. The typical uncertainty associated with the 
data is shown above the extinction vector, which is drawn at the 
bottom right-hand corner. 

the values are compatible with the composite track, would 
typically have a young population of about or younger than 
4 Myr. Were this to be the case, the extinction would be 
negligible and the equivalent widths would be about 200 A. 

- Had we included some Ha contribution from the neigh- 
boring zones then the Ha flux measured would have been 
overestimated. The Ha flux drops considerably as age 
increases across the age interval we consider (by up to more 
than one order of magnitude), but if the Ha flux is weaker 
the mass of the young population decreases accordingly. 
The mass estimated using the red filter can vary by up to 
a factor of four within the age interval 1-10 Myr, according 
to the stellar population models used in this study. The age 
estimate would be older by no more than 1 Myr, if we had 
overestimated the flux by a half. Thus, this correction would 
in our case be negligible. 

Taking into account these sources of uncertainty, the first- 
order estimates of the age were performed for the young 
population inside a complex, assuming a single burst. We 
estimated the age by considering the extinction value that we 
would need to shift the location of each dot (flux quotients) 
along the evolutionary tracks. The estimated ages range from 
2 to 5 Myr (Age phot', see table [3]), the typical age range in 
which an instantaneous burst of st ar formation display s Wolf- 
Ray et (WR) features in its spectra ( iLeitherer et al.ll 1999b . Given 
the spectral range covered in this study, we would expect to 
see the well-known red WR bump, which is the result of the 
blend of the [Ciii 1 /15698 and [Civ] /15808 broad emission lines 
(iKunth&Schildl 11986) . We tried to find these features in the 



spectra of the complexes, but could find no clear evidence 
of these bumps. This bump is much more difficult to detect 
than the other bump (such as the blue bump) characteristic 
of ^VRs because it is always weaker (e.g., iFernandes et al.l 
l2004l: iLopez-Sanchez & EstebanI l2009l) . Sometimes it is not 
even detected a gain in Wolf-Rayet galaxies w here it was 
previously seen (Lopez- S anchez & EstebanI 12009). Therefore, 
the non-detection of WR features does not definitely exclude the 
hypothesis that the young population in our complexes spans the 
age range 2-5 Myr. 

An upper limit to the age can be estimated by using the EW 
(Ha), as it decreases strongly with time. Using the EW peaks 
(see table [21), we estimate that the age is in the range 5-10 Myr 
(Age^^; see table O. In many cases, the EW peaks are not 
strong enough to explain population as young as our former 
age estimates. This is unsurprising because within a complex 
and a spaxel itself the underlying older population contributes 
significantly to the Ha continuum but not to the emission line. 
The contribution of this continuum minimizes the total EW 
within a complex. The broad-band flux of the knots typically 
represents between 1 % and 40% of the total flux within the area 
of the complex. We consider for instance a 3 Myr-old population 
whose broad-band flux represents 10% of its overall broad-band 
flux within the complex. According to the models used in this 
paper, the equivalent width of the young population without 
contamination would be -1000 A. However, the contamination 
of a 1 Gyr-old population would diminish the equivalent width 
to a measured value of 70 A for the whole complex, the same 
value as a single population of 6 Myr without contamination. 
Therefore, with the observed EW we can only set an upper limit 
to the age in each case. 

As a conservative approach, instead of directly using the 
youngest estimates to derive the mass of the young population, 
we use the average of the two estimates, AgCphot and Age^^. 
Thus, once the age and extinction were estimated, the mass 
was directly obtained via the extinction-corrected F814W 
magnitude. 

Under these assumptions, the derived mass of the young pop- 
ulation of the selected Ha complexes is between 10^^ M© and 
10^^ Mq , with the exception of three complexes for which 
the derived mass is about or close to 10^ M0(M[/]; table O. 
These complexes have extinctions of Ay ~ 1-2 mag, with the 
exception of CI in IRAS 14348-1447 (Ay = 4.2 mag). The 
uncertainty in this mass is typically smaller than a factor of 
two. This rather small uncertainty is expected, since during the 
first 1-7 Myr of the starburst the broad-band luminosities do not 
change significantly. However, if we compute the mass using 
the Ha luminosity (M[Ha] ; table [3]) the uncertainties increase 
considerably. This is also expected because the Ha flux for a 
population of a given mass evolves significantly with the age 
of the starburst. The few cases in which both measurements are 
incompatible suggest that the age of the burst is closer to the 
youngest value (AgCphot)- The corresponding M[Ha] would be 
lower toward a value similar to M[/] . 

An old population of a few Gyr (prior to the interaction, 
and in the parent galaxy) in TDG candidates that contributes 
to most of the stellar mass has been reported (e.g. lSheen et all 
,,2009.) . Large Hi reserv oirs have also been found in a few 
TDGs (IDuc et al.l l2007l) that can sustain star formation on a 
Gyr scale. We have seen previously that the colors of some 
complexes can be compatible with a composite population 
where a 1 Gyr-old population is 99% more massive than a young 
burst. Thus, we consider it worth investigating this possibility 
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Fig. 7: Velocity dispersion vs. estimated effective radius. Dots 
correspond to the relation for our selected complexes. The 
typical size for errors is shown in the bottom right corner. The 
lines show the fit for extragalactic Hii regions (continuous), 
elliptical galaxies (da shed) and globular + ellipti cal galaxies 
(dotted) obtained by iTerlevich & Melnickl (Il98lh . The other 
symbols represent diff'erent samples of dynamically hot systems: 
open triangles in blue, intermediate and giant ellipticals 
(Bender et al.l 119921); open sq uares in red, dwarf ellipticals 
(iBender et al . 1992; Geha et al j|20Q3l) : crosses in yellow, dwarf 
spheroidals (Bender et al. .19921): plusse s in green, the TDG 
candidates in Monreal-Ibero et al. (l2QQ7h : diamonds in orange , 
massive globular clusters in NGC 5128 (Martini & Ho" 2004); 
inverse triangles, g lobular clusters in the Gal axy (Trager et al. 
119931) and in M31 (iDubath & GrillmaidlT9^; an d asterisks in 
pink, clusters in the Antennae (iMengel et al 



also here. Table [3] shows the estimates of the age of the young 
(Agecyoung) and the mass of both populations (Mcyoung and 
Mcoid)- The stellar mass of the complexes generally ranges from 
10^^ M© to 10^^ M© , an order of magnitude higher than the 
previous estimates. 

5.2.2. Conditions for self-gravitation 

Measuring velocity gradients helps us assess whether 
there is either any indepen dent rotation (IWeilbacher et al.l 
l2002t iBournaud^etaD l2004 [2008 ) or outflows. However, 
given the spatial resolution achieved with the IFS instruments, 
we cannot in general resolve any velocity field across the 
extranuclear condensations. In this case, other indirect methods 
must be used to assess whether the internal motions of our 
co mplexes can affect their stabi l ity. 

Ilgl esias-Paramo & VilchezI (120011) established an empirical 
luminosity criterion (L(Ha) > 10^^ erg s"^) that must be reached 
by their Hii complexes to ensure self-gravitation. This criterion, 
however, has not been supported by any kinematic study. More 
than 50% of our own complexes fulfill this criterion, as shown in 
table [21 If we corrected for internal extinction, only three would 
be less luminous than 10^^ erg s"^ 

Another method that can be used to establish whether 
the selected complexes are stable is to study their location 
in well-known empirical correlations followed by other self- 
gravitating entities such as elliptical galaxies, bulges of spiral 



Fig. 8: Velocity dispersion vs. Ha luminosity corrected for in- 
ternal extinction and metallicity effects. The lines show different 
fits for extragalactic giant Hii regions which are approximately in 
virial equilibrium, given by: Terlevich & Melnick (1981) (so lid 
line), Relaiio etaP tOOSh (dotted line) and lRozas etal.l (l2006al) 
(dashed line). 



galaxies, globular clusters, and/or giant Hii regions: the radius- 
velocity dispersi on and the luminosity-velocity dispersion 
relations ( Terlevic h & Melnickl Fl 9 8 1 ) . Here, we consider these 
relations because they provide more reliable constraints than the 
Ha luminosity of the complexes. We plot these relations for our 
complexes in Figs. [7] and [H 

In the radius-velocity dispersion diagram, we superimpose 
the data of samples of dynamically hot systems such as massive 
globular clusters, dwarf elliptical and spheroidal galaxies, 
intermediate and giant ellipticals, and TDG candidates. Within 
all of the uncertainties, the velocity dispersions of six of the 
complexes selected here are too high to ensure self-gravitation. 

For a giant Hii region, the lower envelope in the log Lua- 
logcr plane is closely represented by a straight line of the 
form Lua= c -h F x logcr, where cr refers to the velocity 
dispersion, the constant c ranges from 34.5 to 36.8, and F is 
between 2 and 4 (see line fits and references in Fig. [8]). This 
is known as the luminosity-velocity (L-cr) dispersion relation, 
and implies that the giant extragalactic Hii regions above the 
envelope are gravitationally bound complexes of stars and gas 
and that the widths of the nebular emission lines reflect the 
motions of discrete ionized gas clouds in the gravitational field 
of the underlying stellar and gaseous mass ([Terlevich & MelnicS 
[I98l[) . The kinematics of the Hii regions below the envelope 
may be dominated by processes other than gravitati onal motions , 
such as stellar winds and expanding shells (Relaflo et an[2005[: 
[Relafio & Beckmad [20051: [Rozas et al.[[2006a[[b[) . That is. if the 
relation for giant Hii regions lies above the envolope, then 
the non-gravitational processes are identified with broad, low- 
intensity components that do not considerably affect the physical 
properties of the main spectral component (i.e., gravitational 
motions). Thus, although the velocity of these fenomena could 
be similar or even larger than the cr measured in this paper, 
if the relation for our complexes is above the envelope, the 
spectral component corresponding to the gravitational motions 
will probably not be affected. 
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In any case, we estimated the expansion speeds of bubbles 
in the interstellar medium, so as to quantify how broad these 
non-gravitational components are. We evaluated the equivalent 
number of 03 (V) stars using the Ha lumino s ity of a given 
complex and the value given by IVacca et aP (Il996h for the 
ionizing photon output of an 03 (V) star. We then derived the 
ki netic energy inpu t from the stellar winds, using the estimate 
by iLeithereif (1 1998 1) of the wind luminosity for an 03 (V) star, 
integrated over the time the star is on the main sequence. 
Assuming a typical average electron density of 100 cm"^, 
obsevedin (U)LIRGs ( Veilleux et al. 1995, 1999), the expansion 
speeds are expected to be (Lamers & CassineUi 1999) between 
several km s"^ for the least luminous complexes and less than 30 
km s"^ for the most luminous ones. 

Given the metallicity range of our selected complexes, 
we applied a metallicity correction to the Ha luminosity, in 
such a way that Alog (F(Ha))=127z - 1.17 (Terle vich & Melnickl 
11981) , where z denotes the metallicity value (z=0.02 for 
solar metallicity). Within the uncertainties, all the selected 
complexes but two are consistent with the line fits or above the 
envelope (Fig. [8]). Interestingly, these two outlying complexes 
are consistent with the radius- velocity dispersion correlation. 
Although the values for another complex are below the envelope, 
we note that its Ha luminosity was not corrected because its 
metallicity is unknown. Thus, no complex fails both criteria at 
the same time. According to these results, self-gravitation should 
be ensured for more than a half of the complexes for which a 
determination of the velocity dispersion was achieved. For the 
remainder, self-gravitation is neither guaranteed nor excluded 
because at least they generally satisfy one of the two criteria 
used. 



5.2.3. Dynamical mass estimates 

Upper limits to the dynamical masses (Mdyn) of the selected 
complexes were derived under the following assumptions: (i) the 
systems are spherically symmetric; (ii) they are gravitationally 
bound; and (iii) they have isotropic velocity distributions 
[(T^(total)=3X(T^^^], where crios is the line-of-sight velocity 
dispersion. In the previous section, we have seen that more than 
a half of these complexes are likely to be gravitationally bound. 

The mass of a virialized stellar system is given by two 
parameters, its velocity dispersion (cr) and its half mass radius 
(vhm)' The half-mass radius is not an observable and cannot be 
measured directly. It has to be inferred indirectly by measuring 
the half-light radius (i.e., r^ff). For a cluster-like object this mass 
is given by 



Mdyn/Mo = rj 



G 



(1) 



where the reff is given in pc, ctlos in km s \ the g ravitational 
constant 0=4.3x10"^ pc M© (kms-i)^ dSpitzerlll987b . and 
T] = 9.75 for a wide range of light profiles. However, some 
studies have shown that rj is not a constant, and can vary 
with time (77 ~ 3 - 10), depending on, for example, the degree 
of mass se greg ation and the binary fraction of the cluster 
(iFleck et al.ll20Q6l: lKQuwenhoven & de Griisll2008l) . Since we do 
not have this information, we use 77 = 9.75. 

Since the star-forming knots ha ve a cluster-like light profile 
(see 'Miralles-Ca ballero et al ] 1201 lb the use of Eq. □ to derive 
the dynamics of the selected complexes seems appropriate. The 
corresponding masses given in table O are between a few 10^ 
and a few 10^ M©, with an uncertainty of a factor of 2-3. 



5.2.4. Total mass: discrepancy between the stellar and the 
dynamical masses 

It is unclear where to draw the line between a high-mass 
young super cluster and a low-mass dwarf galaxy in terms 
of stellar mass. We derive masses only for young population, 
and those complexes less massive than 10^ M© are comparable 
to the brightest young clusters or complexes hosting recent 
star formati on in other less luminou s interacting galaxies 
(e.g., Mengel et all l2008l; iBa stian et a lJ ^2009). However, does 
the estimated mass of the young stellar population represent 
the total mass of the candidates? Only a few studies have 
reported stellar masses for TDGs. The dynamical mass of our 
complexes is in the range of the typical total baryonic mass of 
a TDG (10^-10^ Mc, ; e.g.. Due et al. 2004; Higdon etaD[2QQ6l: 
iHancock et al.ll2009l) . However, the ratio of the dynamical to the 
most massive photometric mass in our complexes is suspiciously 
high (Mdyn/m[M;] ~ 50-1000). Several factors can explain this 
ratio. 

If we assume that the stellar mass represents the total mass 
of a given complex, this mass might be underestimated, because: 

- Although the most recent bursts of stars are responsible 
for the bulk of the ionization of the gas - hence of the 
emitted flux- the kinematics of the ionizing gas is afl'ected 
by those of previous stellar generations. In regions of mixed 
populations, the light is generally dominated by the young 
stars a nd the mass by th e older population. As outlined 
before, ISheen et aP (l2009l) derived a total stellar mass of 
3.1x10^ M© for one of their TDG candidates, a fraction of 
no more than 2% being contained in populations younger 
than 6 Myr. The mass-to-light (M/L) of young (e.g., 5 
Myr) and old (e.g., 1 Gyr) populations difl'ers by about a 
factor of 100. The Ha to the broad-band luminosity colors 
of some of our candidates are also consistent with those 
of composite populations. Our complexes might indeed 
consist of young star-forming knots and evolved population 
from the complexes themselves and the parent galaxy that 
has been accreted by the complex. We estimated the total 
stellar mass of the complex if practically all the F814W 
flux measured inside the whole complex (minus the flux 
from the knots) belonged to it and came from a 1 Gyr-old 
population. This total stellar mass approaches significantly 
to the dynamical mass measured for most of the complexes. 

- The total stellar population does not normally represent 
the total baryonic mass (gas -\- stellar), and we do not have 
access to the mass of the gas. Normally, the efficiency of star 
formation in molecular clouds is very low (typically about 
1%). In (U)LIRGs the efficiency can be higher, but the mass 
of a molecular cloud undergoing a starburst episode may still 
be higher than the mass of the starb urst by more than one 
order of magnitude. The candidates in lHancock et al.l (l2009h 
have total stellar masses of 1-7x10^ M©, the HI mass being 
6x10^ M©. Studies based on radio Hi observations have 
found that the gas mass of a TDG with a typical diameter 
of a few kpcs lies wit hin the range between a few 10^ and 
several 10^ M© (e.g., lDuc et al.ll2000l: iBournaud etaDl2004l: 
iDuc et al]l2007h . 



On the other hand, if we assume that the dynamical mass 
represents the total mass of a given complex, we might be 
overestimating it, because: 
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- In each case, Mdyn represents normally an upper limit to the 
dynamical mass. For instance, given the complex dynamics 
in interacting systems, high velocity dispersions might not 
necessarily indicate high masses, but might alternatively 
represent tidal-flows induced from the merging process and 
strong winds from LINER -like regions (IColina et al.ll2005l: 
iMonreal-Ibero et a02QlC)l) . However, at least for the systems 
that we previously observed with VIMOS, these strong 
winds may not afl'ect sign ificantly the de termination of the 
velocity dispersions (Monreal -Ibero et ani2QlQ ). 

- Dark matter might contribute significantly to the total TDG 
mass. However, TDGs are unlikely to contain a large amount 
of dark matter (Barnes & Hernquist 1992) because their 
material is drawn from the spiral disk, while the dark matter 
is thought to surround the galaxy in an extended halo. 

On the basis of all these arguments, it is clear that the 
total masses of these complexes is somewhere between the 
photometric and the dynamical masses. Other independent 
criteria also show that more than a half of the selected complexes 
are in virial equilibrium, thus the velocity dispersion traces 
more li kely th e dynamical mass rather than other processes (see 
section [522I)- Given all the caveats and uncertainties previously 
considered when determining the dynamical and stellar masses, 
we assume as a mass criterion that the complexes with 
photometric mass estimates using the broad-band luminosities 
(single or composite population) that are compatible with 
10^ M© are likely to have a sufficient total mass to become a 
tidal dwarf galaxy. This assumption should obviously be verified 
using multi-wavelength photometric and spectroscopic data. 
Only six of the selected complexes (plus the three candidates in 
IRAS 16007-1-3743) fulfill this mass criterion. Had we assumed 
that a complex with a dynamical mass compatible with 10^ M© 
or higher has sufficient mass to become a TDG, then 15 of the 
selected complexes would fulfill the mass criterion. 

5.3. Are the selected Ha complexes unaffected by the forces 
from the parent galaxy? 

5.3.1. Tidal forces 

If the complex is massive enough to constitute a dwarf galaxy 
and gravitationally bound, its fate basica lly depends on the 
ratio of its mass to the so-called tidal mass (iBinney & Tremaing 
119871: lMende~ de Oliveira et al ]l2Q0r). The tidal mass condition 
will tell us whether it is massive enough to survive the tidal 
forces exerted by the parent galaxy. The tidal mass is defined 
as 

Mtid = 3xM|^j , (2) 

where M refers to the mass of the parent galaxy, R is the radius 
of the complex (here the size of the most luminous knot in a 
complex, as derived in MCll), and D is the distance to the 
parent galaxy. This equation is valid when the size (R) of a 
certain region is small in comparison with the distance (D) to 
the parent galaxy, which is the case for all of the complexes 
(the ratio "size/projected distance" is typically below 0.04). If 
the tidal mass of an object is lower than the total mass, then the 
object is unaffected by the forces applied by the parent galaxy. 

In general, the gravitational potential of interacting galaxies, 
as in ULIRGs, is a complex function of the mass distribution 
of the system, which evolves with time. Nevertheless, as a first 
approximation, we assumed that the distribution is dominated 



by the masses of the main bodies of the system. Two different 
approaches were applied in this study: the gravitational potential 
depends either on: i) the nearest galaxy (M||^^0 or ii) a point mass 
in the mass center and the total mass of the system (M™). The 
highest value obtained using the different approaches was taken 
to compare with the dynamical mass of the complex. 

In each case, we used the measured projected distances Jnear 
and dcM, respectively (see table [T]). Since the projected distance 
is always smaller than the real one, we actually estimated upper 
limits to the real tidal masses. To derive the tidal mass, we need 
to know the mass of the parent galaxy, which in each case is 
identified with its dynamical mass. Under the same hypothesis 
as in section [5^31 of a virialized system, the dynamical mass of 
a galaxy-like object can be obtained as 

Mdyn/Mo = mlO^rhm X , (3) 

where, the half-mass radius, Vhm, is given in kpc, cr in km s"^ 
and the factor m is a dimensionless function of the assumed 
mass distribution and ranges from 1.4 for a King stellar mass 
distribution that adequ ately represents ellipticals (jBender et all 
ll992l:Pl^coni et al1l200Z) to 1.75 for a poly tropic sphere with a 
density index covering a range of values ( Spitzer 1987) a nd 2.0 9 
for a de Vaucoluleurs mass distribution (Combes et al1 ll995h . 
We assumed that m=1.75, as in Colinaetal. (2005). In each 
case, we used the average value of the velocity dispersion at the 
peak and nearby surrounding spaxels, covering about a radius of 
V\ Finally, the effective radii (from which we infer Yhm) of the 
parent galaxies were derived using the //-band NICMOS and 
WFC3 images of the galaxies and fitti ng thei r two-dimensional 
fight profiles with GALFIT (IPeng et al.ll2QQl iPeng et al.ll2010l) . 
In many cases, the NICMOS images do not cover the whole 
extent of the galaxies, which meant that we had to use the values 
given in Arribas et al. (in prep.). 

In summary, two values (M||^^^ and M|r^) for the tidal 
mass were obtained using Eq. O and the maximum value was 
compared to the dynamical mass of the selected complexes 
derived in section \5?23[ The ratio Mdyn/Mtid is in general > 10 
(see tableO, which ensures that the complexes are unaffected by 
the tidal forces exerted by the parent galaxy. Only one complex 
fails this condition and another one is close to failing because of 
large uncertainties in its dynamical mass estimate. 

5.3.2. Escape velocity 

If the complex does not have a large enough relative velocity 
(Vrei) with respect to the galaxy to escape, it might still fall back 
towards the center of the system because of the gravitational 
force exerted by the parent galaxy. It is interesting to consider 
whether a complex does indeed exceed the effective escape 
velocity (Vesc)- This criterion, however, should not have the same 
weight as the others since it is subject to many uncertainties: 
(i) the real distance of the complex is always larger than 
the projected one; (ii) only one component of the velocity is 
measured and there is no information about the movements in 
the plane of the sky; and (iii) two possibilities for the relative 
movements between the region and the system are always 
possible, since for a given configuration it is impossible to 
determine whether the complex is closer or further away from 
the observer than the mass center. We include this criterion in 
this study for completeness. 

For simplicity, we assume that the gravitational potential 
is created by a point mass representing the total mass of the 
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system, that is located at the mass center. The ratio |Vescl/Vrei 
(last column in table [3]) is smaller than 1 for 7 of the selected 
complexes out of the 22, that is, they could escape. If we take 
into account projection effects statistically, a complex escapes if 
Ivreil - Ivesc X cos(7i/4)\ > 0. In this case, an additional complex 
satisfies the criterion. 



5.4. How common is TDG formation in (U)LIRGs? 
5.4.1. TDG candidates in (U)LIRGs 

None of the developed criteria can help us determine wheter 
a given complex will survive the merging process. Projection 
eff'ects and observational constraints, especially the need for 
higher angular resolution spectroscopic data (which would allow 
us to search for velocity gradients), ensure that it is difficult to 
assess the fate of the selected complexes. In any case, we can 
investigate which have the higher probabilities of surviving as a 
TDG, based on the fulfillment of a few or most of the criteria 
viewed in this study. The results are summarized in table IH 

We assigned different weights to the diverse criteria we 
studied and derived a probability that a certain complex could 
survive as a TDG by adding the weights for the criteria that 
it achieves. The mass criterion is considered to be the most 
important since, even if it is self-gravitating, if it does not have 
enough mass the complex could be either a bound super cluster 
or any other entity rather than a low-mass dwarf galaxy. For 
this reason, it is assigned a weight of 30%. The criteria with 
the least importance (each having a weight of 10%) correspond 
to the Ha luminosity and the escape velocity ones, since the 
former comes from empirical considerations and the latter is the 
least reliable of all. We assigned a weight of 20% to the Mtid 
versus (vs.) Mdyn criterion. Finally, we considered a somewhat 
lower weight (15%) for each of the self-gravitating criterion (cr 
vs. Teff and cr vs. Lua), since they prove the same condition. In 
practice, the complexes that satisfy most of the criteria have the 
highest probabilities of being TDG progenitors. 

Once the percentage was computed based on which criteria 
the complexes fulfill, different probabilities that a complex 
constitutes a TDGs were defined: low (prob < 40%), medium 
(40 < prob < 60), medium-high (60 < prob < 80), and high (prob 
> 80). Only three complexes (counting IRAS 16007 -h3743 have 
a high probability. A total of 6 of the 22 selected complexes (9 of 
25 if we include those from IRAS 16007-h3743) have medium- 
high or high probabilities of being TDG progenitors, and from 
now on we consider these as our TDG candidates. 

We detect candidates in LIRGs and ULIRGs, although the 
three candidates with a high probability are only found in 
ULIRGs. This suggests that TDG production may be more 
efficient in systems with higher infrared luminosities. However, 
this statistic is of very limited robustness because of the low 
number of candidates found. If we consider only systems for 
which IFS data are available and located at a distance (i.e., Z)l > 
130 Mpc) where at least a galactocentric radius of 10-15 kpc 
is covered with the spectroscopic data, we obtain a production 
rate of about 0.3 candidates per system for the (U)LIRG 
clas s. Other TDGs have been observed at la rger distances (>30 
kpc: ISheen et al ] |2009I: iHancock et ani2009h not covered by the 
FoV of our IFS data. The detection of bright and blue knots (i.e., 
Mp^uw > -12.5 and MF/^2>5w-MFum < 0-5) in the tidal tails and 
their tips in the ACS images of some (U)LIRGs (regions not 
cover ed with the IFS data; see Fig. 8 in Miralles-Caballero et al.l 
l201lh . indicative of a young population more massive than 



10^ Mq, may increase the number of candidates for the (U)LIRG 
class. 

5.4.2. Dynamical evolution of the TDG candidates 

The number of candidates in systems undergoing the early 
phases of the interaction process (i.e., phases I-II and II) is 
7, significantly larger than the 2 detected in more evolved 
systems. Normalized to the total number of systems in early 
and advanced phases, we obtain 0.5 and 0.13 candidates per 
system, respectively. Therefore, TDGs are more likely to be 
formed during the first phases of the interaction in (U)LIRGs. 

To analyze the meaning of this trend we need to consider 
both the life expectancy of these candidates as well as their 
detectability during their lifetime. Our result does not mean 
that the total number of candidates in each phase decreases 
at all. Although large Hi reservoirs have been found in a few 
TDGs (I Duc et al] 12007 ) that would be able of sustaining star 
formation on a Gyr scale, it is unknown whether these are 
common. Our candidates are indeed Hof -selected. We might 
have been unable to identify candidate objects as old as 20 Myr 
and older, for which their Ha emission is undetectable. There 
is a score of knots in the initial sample of 32 galaxies at a 
projected galactocentric distance larger than 2.5 kpc, with colors 
(Mf435w-Mfu4w > 1-5) and luminosities (Mfu4w < -15), such 
that if the population were about 100-1000 Myr old the stellar 
mass would be higher than lO'^-lO^ M©, close to or similar to the 
tot al stellar mass of o b served TDG candidates. 

iBournaud & Dud (l2006h ivestigated a set of 96 N-body 
simulations of colliding galaxies with various mass ratios 
and encounter geometries, including gas dynamics and star 
formation. They investigated the dynamical evolution of the 
TDG candidates found in the various simulations up to 
t = 2 Gyr after the first pericenter of the relative orbit of the 
two galaxies. On the basis of the compariso n with a higher 
resolution simulation of a major merger by iBournaud et al.l 
(2008), we can roughly assign different dynamical times after 
the first pericenter for the d ifferent phases of interaction defined 
in Miralles-Caballe ro" et al .^ (2011). In this way, the first phases 
of the interaction (I-II and III) are likely to occurr during the 
first 400 Myr of the interaction process. If we consider only 
the observed candidates in (U)LIRGs during these early phases 
and these systems for which we have available IFS data that 
cover at least 10-15 kpc of galactocentric radius we derive a 
production rat e of about 0.6 candidate s per system. Although the 
simulations of IBournaud & Dud ([2006) continue up to t = 2 Gyr, 
they fit the number of their TDG candidates versus time with 
an exponential decay and a lifetime of 2.5 Gyr. They estimate 
the number of candidates that survive (long-lived candidates 
that neither fall back nor lose a large fraction of their mass) 
after 10 Gyr corresponds to 20% of the TDG candidates formed 
during the early phases of the interaction (i.e., t < 500 Myr 
after the first pericenter). According to this percentage, we 
estimate a production rate of about 0.1 (20% of 0.6) long-lived 
TDG candidates per system for the (U)LIRG class. 

The average production rate of 0.1 long-lived TDGs per 
system can either become lower or higher depending on 
follow-up studies with IFS capable of providing higher spatial 
resolution (i.e., adaptive optics assisted systems) and covering 
a FoV that will even allow us to study galactocentric distances 
greater than 30 kpc. Modeling and observations claim that 
the most prominent TDGs are found along the tidal tails 
and in particular at t h eir tips (e.g., iHi bbard & Mihos 1995; 
iDuc & Mirabeil 119981: IWeilbacher etaLI >2000: .Higdon et all 
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120061: iDuc et alJ 120071: iBournaud et aP 120081: iHancock et aP 
'2009'). Thus, it would be unsurprising that when covering the 
whole field of the local (U)LIRG systems with IFS facilities, 
more TDG candidates will be found. In addition, kinematic data 
of higher spatial resolution would allow us to determine more 
reliably whether the candidates are bound. 

5.4.3. Implications for TDG formation at high-z 

After estimating the production rate of long-lived TDG candi- 
dates per system for the (U)LIRG class, we can now roughly 
estimate whether satellites of tidal origin are common, and more 
specifically, the contribution of TDGs in the early Universe 
to the dwarf population we see in the local universe. Local 
(U)LIRGs represent an appropriate class of objects to study 
the importance of TDGs at high redshift because: (i) they 
are maj or contributors to the star formation rate density at 
z ~ 1-2 (iPerez-Gonzalez et al.|[2005 ): and (ii) ULIRGs resemble 
the sub-millimeter galaxies detected at higher redshifts in the 
sense that they are merging systems with e xtremely high rates 
of star formatio n (Chapman et ^ 120031: iFraver et al.1 120031: 
lEngel et al.ll2010h. ^ 

^Qk azaki & Tanig uchil (|2000l) studied a scenario in which 
galaxy interactions and/or merger events act as the dominant 
formation mechanism of dwarf galaxies in any environment 
and these interactions occur in the context of the hierarchical 
structure formation in the Universe. They claimed that all dwarfs 
in the Universe might have a tidal origin. On the basis of this 
scenario, their statement assumed a production rate of 1 to 2 
tidal dwarfs per merger, which has a lifetime of at least 10 Gyr. 
The production rate we estimate for long-lived candidates in 
(U)LIRGs is a factor of 10-20 lower, such that their contribution 
to the overall dwarf population is only 10-5%. Therefore, we 
do not find strong evidence that all the local dwarf galaxy 
population have a tidal origin. This result is consistent with 
recent estimates implying that the contribu tion of TDGs to 
the overall dwa rf population is insignific ant (iBournaud & Dud 
IIOOS |^nitani20Tl'; 

The study by Okazaki & Tani guchi1 (l2000l) covers a wide 
range of mass ratios for the interacting systems. In our sample, 
the systems with more than one nucleus span mass ratios from 
1:1.7 to 1:3.6, which is the most favorabl e range for long- 
lived T DGs, according to the simulations of Bournaud & Dud 
(I2006h . Were high-z interac tions to have similar mass rat ios, the 
production rate estimated by lOkazaki & Taniguchil (l2000l) would 
be higher than the predicted 1-2 tidal dwarf per merger. Thus, 
our production rate would be a factor of more than 10-20 higher, 
and the contribution of TDGs to the overall d warf population 
would become even lower than 10-5%. Rodi ghiero et al.l (1201 ll) 
claimed that mergers at z ~ 2 have a very low contribution 
(~ 10%) to the cosmic star formation rate density. Were this 
result to be confirmed, the contribution of TDGs to the overall 
dwarf population would still be lower, to a negligible value. 

6. Conclusions 

We have combined high angular resolution HST images 
with spectroscopic data from IFS facilities to characterize 
Hof-emitting clumps in an initial sample of 27 (U)LIRGs. 
Ou r study has extend e d the search of TDG candidates 
by iMonreal-Ibero et al.l (l2007h by considering both LIRG 
systems and additional and more sensitive images. In particular, 
we have detected a total of 31 extranuclear star-forming 
complexes in 1 1 (U)LIRGs and characterized the main physical 



and kinematic properties of these complexes. Using structural, 
physical, and kinematic information, we have also estimated 
the stellar and dynamical mass content of the complexes and 
studied their likelihood of resisting the eff'ects of internal as well 
as external forces based on diff'erent dynamical tracers. With 
all these parameters, we have identified which complexes have 
the highest probabilities of being long-lived TDG candidates in 
(U)LIRGs. We draw the following conclusions: 

- Located at an average projected distance of 9.3 kpc from 
the nucleus, within the range 3.3-13.4 kpc, the structures 
of the complexes are generally simple, consisting of one 
or a few compact star-forming regions (knots) in the 
HST images, though a few of them reside in a richer cluster 
environment. The complexes have typical B luminosities 
(Mf435w < -10.65), sizes (from few hundreds of pc to about 
2 kpc), and B - I colors (Mf435w-^FSi4w^ 1-0) that are 
similar to those observed in giant Hii regions and dwarf-like 
objects. The relatively high metallicities derived, of normally 
Zq-Zq/3 (independent on the luminosity of the complex), 
reflect the mixing of the metal content in interacting 
environments, as observed in a few extragalactic Hii regions 
and TDG candidates. 

- The measured Ha luminosities of the complexes are compa- 
rable to those of extremely luminous Ha complexes in nearby 
systems, giant extragalactic Hii regions, TDG candidates, and 
normal dwarfs. Their sizes seem to be more related to the 
nature of the system. In our case, many complexes have sizes 
typical of dwarf-like objects. 

- Twenty-two complexes with very simple structures were 
selected to study their nature as TDGs. The stellar masses 
derived for these complexes using the Ha and broad-band 
luminosities and equivalent widths of a single burst or a 
composite population (young burst -h 1 Gyr) range from a few 
10"^ to 10^ Mq. Nevertheless, it is typically below 10^-^ M©, 
which is the lower limit to the total mass. In contrast, the 
complexes have dynamical masses that are a factor of 50-1000 
higher, which provides an upper limit to the total mass. 

- A total of 9 complexes, namely TDG candidates, have the 
highest probabilities of becoming dwarf galaxies, at least up 
to 1-2 Gyr. They are probably massive enough and satisfy 
most of our criteria for self-gravitation (i.e., the position in the 
radius-cr and the luminosity-cr planes) and resistance to the 
forces from the parent galaxy. We have found evidence that 
their formation takes place more often in early phases of the 
interaction. 

- When we consider only systems for which the IFS data cover a 
significant fraction of the whole system, the production rate of 
candidates averages about 0.3 per system. This rate is expected 
to decrease to 0.1 for a long-lived 10 Gyr dwarf, according to 
recent galaxy merger simulations. Were this to be the case, 
fewer than 5-10% of the general dwarf satellite population 
could be of tidal origin. 
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Table 1: Photometric properties of the identified star-forming complexes: distances, magnitudes, and sizes 
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Notes. Col (1): object designation in the IRAS Point and Faint Source catalogues. In case of several pointings taken with VIMOS 
data, the orientation (N: North, S: South, SE: South-East, C: Center) of the centered nucleus in that pointing is given. Col(2): Number 
of nuclei in the system (1-3) and interaction phase (IP) that the system is undergoing according to the classification scheme defined 
in Mir alles-Caballero et al.l (|201 1), where I-II indicates first approach. III pre-merger and IV merger phases, respectively. Col (3): 
identified complex. Col (4): number of knots per complex. Col (5): projected distance to the nearest galaxy. Col (6): projected 
distance to the mass center of the system. Col (7): Absolute magnitude, M/7435, not corrected for internal extinction. Col (8): derived 
photometric color of the complex. Col (9): r^^ of the brightest knot inside the complex. Col (10): equivalent radius, derived using 
the size of the knots (see text). Col (11): radius of the complex measured on the Ha map. 
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Table 2: Spectral observables, metallicity, and dynamical parameters of the identified star-forming complexes 
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1 


7.0 


6.0 


56.6 


242 


4.3 


46 


196 


8.53 


08572+3915 SE 


1 


15.8 


13.6 


52.2 


183 


3.5 


74 


148 


8.54 




2 


5.8 


5.0 


4.1 


85 


21.0 


62 


207 


8.50 




3 


1.0 


0.8 










228 


8.53 




4 


1.0 


0.8 


8.8 


28 


3.2 






8.56 


F10038-3338 


1 


12.6 


4.1 


86.0 


374 


4.4 


19 


19 


8.68 




2 


0.6 


0.2 


34.3 


74 


2.2 




-11 


8.32 




3 


2.0 


0.6 








17 


-32 






4 


2.2 


0.7 








27 


4 


8.24 


12112+0305 


1 


44.8 


60.6 


71.1 


252 


3.6 


72 


96 


8.59 




2 


11.9 


16.1 


26.3 


41 


1.6 


81 


-204 


8.81 




3 


8.7 


11.8 


11.5 


36 


3.1 




-349 


8.59 




4 


4.9 


6.7 


6.9 


53 


7.7 




-266 


8.62 


14348-1447 


1 


29.5 


65.7 


88.9 


149 


1.7 


54 


106 


8.63 


15250+3609 


1 


8.1 


5.7 


168.9 


258 


1.5 


69 


124 


8.61 


F18093-5744 N 


1 


32.8 


2.6 


48.1 


213 


4.4 


23 


-76 


8.75 




2 


71.8 


5.6 


67.0 


263 


3.9 


24 


60 


8.77 




3 


9.5 


0.7 


15.8 


21 


1.4 


39 


48 


8.90 


Fl 8093-5744 C 


1 


3.8 


0.3 


26.0 


57 


2.2 


18 


119 


8.62 


23128-5919 


1 


5.8 


2.8 


32.0 


51 


1.6 


69 


2 




16007+3743 


Rl 


12 


57.1 


40 






61 


-41 


8.7 




R2 


35.5 


189.5 


234 






76 


92 


8.7 




R3 


5.85 


8.8 


80 






85 


338 


8.7 



N otes. Col (1): object designation as in table[T] Col (2): identified complex. Col (3): Observed Ha flux. As stated in lGarcia-Marm et a n (l2009bh and 
in "Rodriguez Zaur m et al.l (1201 Ih . typical uncertainties in the absolute flux calibration are between 15% and 20%. Col (4): derived Ha luminosity 
with no correction for internal extinction, except for IRAS 16007+3743. Col (5): computed equivalent width (EW) from the Ha and EW maps 
(see text). Col (6): Spaxel with the largest value of the EW. Col (7): ratio of the peak of the EW to the derived EW for the whole complex. Col (8): 
central velocity dispersion. Col (9): relative velocity with respect to the mass center of the system. Col (10): adopted metallicity for the complex, 
derived using two different line-ratio calibrators (see text). 



Table 3: Derived characteristics of the stellar populations and dynamics of the identified star-forming complexes 



Tr> A c 


Complex 
Number 


Agephot 
(Myr) 


AgeEw 
(Myr) 


Ay 

(mags) 






''^g^Cyoung 

(Myr) 


Ave 
(mags) 


A /f 

Mcoid 


-M^Cyoung 


Mdyn 


A 4^near 
Hid 


A/rCM 

Hid 


Mdyn/Mtidal 

(km s-i) 


^esc 


|Vescl/Vrel 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


04315-0840 


1 


3.4 


6.0 


1.7 + 0.3 


4.9 + 0.1 


5.4 ± 0.3 


2.1 


0.2 + 0.4 


5.8 + 0.1 


3.8 + 0.1 


7.5 + 0.3 


6.6 + 0.3 




8.1 


224 


8.6 




2 


4.0 


10.0 


1.5 + 0.4 


4.2 + 0.2 


5.8 + 0.5 


2.8 


0.7 + 0.5 


5.5 + 0.2 


3.5 + 0.2 


7.8 + 0.3 


6.0 + 0.3 




> 10 


235 


2.3 


06076-2139 


1 


2.9 


5.4 


1.1 +0.3 


5.4 + 0.2 


5.9 + 0.2 










8.0 + 0.5 


5.8 + 0.2 




> 10 


96 


0.2 




2 


3.0 


6.0 


1.0 + 0.4 


4.9 + 0.2 


5.4 + 0.3 










8.0 + 0.3 


5.7 + 0.2 




> 10 


96 


0.2 




3 


4.4 


6.4 


0.6 + 0.3 


4.6 + 0.1 


5.1 +0.2 


2.9 


0.0 + 0.0 


5.6 + 0.1 


3.6 + 0.1 


7.9 + 0.3 


4.5 + 0.2 




> 10 


92 


0.2 




4 


3.0 


6.5 


0.8 + 0.4 


4.5 + 0.2 


5.0 + 0.5 










8.0 + 0.3 


4.5 + 0.2 




> 10 


93 


0.2 




5 


2.9 


6.8 


0.4 + 0.2 


4.4 + 0.2 


5.1 +0.4 










7.6 + 0.4 


4.3 + 0.2 




> 10 


97 


0.2 




6 


2.8 


7.0 


1.0 + 0.3 


5.2 + 0.1 


5.9 + 0.8 












5.0 + 0.2 






103 


0.2 


07027-6011 S 


1 
2 


2.1 
3.2 


5.0 


1.2 + 0.3 
0.6 + 0.3 


4.9 + 0.2 
4.5 + 0.1 


5.6 + 0.3 
4.6 + 0.2 










7.9 + 0.4 


6.1 +0.3 
5.6 + 0.3 




> 10 


180 
125 


7.5 
1.7 


08572+3915 N 


1 


2.8 


5.2 


0.7 + 0.3 


5.4 + 0.2 


5.9 + 0.2 










8.7 + 0.3 


7.2 + 0.4 


7.4 + 0.3 


> 10 


233 


1.2 


08572+3915 SE 


3 


4.9 




0.9 + 0.3 


5.1 +0.1 


5.1 +0.2 


3.1 


0.0 + 0.6 


6.1 +0.2 


4.1 +0.2 




6.6 + 0.4 


6.6 + 0.3 




225 


0.99 




4 


4.9 


6.9 


1.4 + 0.4 


5.5 + 0.2 


6.1 +0.2 


3.2 


0.6 + 0.5 


6.5 + 0.2 


4.5 + 0.2 




6.3 + 0.4 


6.5 + 0.3 




207 




F10038-3338 


3 


3.6 




0.7 + 0.3 


4.4 + 0.1 


4.5 + 0.1 


2.7 


0.0 + 0.0 


5.7 + 0.1 


3.7 + 0.1 


7.6 + 0.3 


5.8 + 0.3 




> 10 


117 


3.6 




4 


4.9 




0.9 + 0.3 


5.0 + 0.1 


5.2 + 0.1 


3.0 


0.0 + 0.6 


6.0 + 0.2 


4.0 + 0.2 


8.2 + 0.3 


5.6 + 0.3 




> 10 


113 


26.2 


12112+0305 


1 


3.4 


5.2 


2.2 + 0.3 


7.1 +0.2 


7.3 + 0.2 


2.2 


0.8 + 0.4 


8.0 + 0.1 


6.0 + 0.1 


9.4 + 0.2 


9.2 + 0.3 


9.6 + 0.2 


0.7 


496 


5.1 




4 


2.1 


6.3 


1.1 +0.3 


5.4 + 0.2 


6.1 + 0.3 












6.2 + 0.3 


6.7 + 0.2 




275 


1.03 


14348-1447 


1 


3.3 


5.5 


4.2 + 0.3 


7.7 + 0.1 


8.1 +0.2 


1.5 


2.5 + 0.3 


8.5 + 0.1 


6.5 + 0.1 


9.3 + 0.2 


9.1 + 1.1 


9.2 + 0.3 


1.3 


575 


5.4 


15250+3609 


1 


5.4 


5.2 


2.4 + 0.3 


6.6 + 0.1 


6.7 + 0.2 


3.9 


1.1 +0.4 


7.5 + 0.1 


5.5 + 0.1 


9.3 + 0.3 


7.7 + 0.2 




> 10 


214 


1.7 


Fl 8093-5744 N 


3 


2.9 


7.2 


1.9 + 0.3 


5.0 + 0.2 


5.8 + 0.5 










7.9 + 0.4 


4.1 +0.1 


5.1 +0.2 


> 10 


220 


4.6 


Fl 8093-5744 C 


1 


3.5 


6.2 


1.6 + 0.3 


4.8 + 0.1 


5.2 + 0.3 


2.4 


0.4 + 0.4 


5.7 + 0.1 


3.7 + 0.1 


7.6 + 0.5 


3.6 + 0.3 


5.6 + 0.2 


> 10 


175 


1.5 


23128-5919 


1 


3.3 


6.3 


1.8 + 0.3 


5.7 + 0.2 


6.1+0.5 


2.0 


0.3 + 0.4 


6.5 + 0.2 


4.5 + 0.2 


9.0 + 0.2 


7.1 +0.2 


7.0 + 0.3 


> 10 


222 


85.7 


16007+3743 


Rl 




7.1 


1.7 


8.8 












9.8 


8.2 


8 


> 10 


302 


7.4 




R2 




5.4 


1.2 


8.8 












10 


9.2 


10.1 


0.9 


408 


4.4 




R3 




6.4 


2.3 


7.8 












10.1 


10 


9.1 


1.2 


519 


1.5 



Notes. Col (1): object designation as in table [T] Col (2): identified complex. Col (3): age of the young population, derived using the photometric information (Fig.|6]). Col (4): age of the young 
population using the EW (Hof). Col (5): internal extinction. Col (6): derived mass of the young population using the I magnitude. All masses in the table are given in log (M (M©)). Col (7): derived 
mass of the young population using the Hof flux. Col (8): derived age of the young component, assuming that the candidate consists of a composite population. Col(9): derived internal extinction of 
the composite population. Col (10): derived mass of the old population (i.e., IGyr), assuming a composite population (see text). Col (11): derived mass of the estimated young population, assuming 
a composite population. Col (12): dynamical mass. Col (13): tidal mass assuming the potential is created by the nearest galaxy. Col (14): tidal mass assuming the potential is created by a point mass 
at the mass center of the system. Col (15): ratio of the dynamical to the tidal mass. Col (16): escape velocity. Col (17): ratio of the escape to the relative velocity. 



D. Miralles-Caballero et al.: Extranuclear YLa complexes in low-z (U)LIRGs: Precusors of TDGs? 
Table 4: Summary of the different criteria used to investigate the nature of the complexes 



IRAS Complex Mass L (Ra) cr vs. rgff cr vs. Mtid vs. Mdyn Vesc vs. Vrei Prob 



Number 


04315-0840 


1 


N 


Y 


Y 


Y 


Y 


N 


IVTpHiiim 




2 


N 


N 


N 


Y 


Y 


N 


Low 


06076-2139 


1 


N 


Y 


Y 


Y 


Y 


Y 


IVTpHinm-Vii oVx 




2 


N 


Y 


Y 


Y 


Y 


Y 


TVTprliiitTi-lii ctIi 




3 


N 


N 


Y 


Y 


Y 


Y 


IVTprliiirn 




4 


N 


N 


Y 


N 


Y 


Y 


A/TpHiiim 




5 


N 


Y 


Y 


Y 


Y 


Y 


A/TpHiiim -hi o\\ 




5 


7 


Y 


7 


7 


7 


Y 


Low 


07097-601 1 S 

\J 1 v/Z/ / \J\J 1 1 O 




N 


Y 


N 


Y 


Y 


N 


IVTpHiiiTTi 

1V1CU.1 Ulil 




2 


N 


Y 


7 


7 


7 


N 


Low 


08572+3915 N 


1 


N 


Y 


Y 


Y 


Y 


N 


Medium 


08572+3915 SE 


3 


? 


Y 


7 


7 


Y 


Y 


Low 




4 


Y 


Y 


7 


7 


7 


7 


Low 


F10038-3338 


3 


N 


Y 


Y 


7 


Y 


N 


Medium 




4 


Y 


Y 


Y 


N 


Y 


N 


Medium-high 


12112+0305 


1 


Y 


Y 


N 


Y 


N 


N 


Medium 




4 


? 


Y 


7 


7 


7 


Y 


Low 


14348-1447 


1 


Y 


Y 


Y 


Y 


Y 


N 


High 


15250+3609 


1 


Y 


Y 


N 


Y 


Y 


N 


Medium-high 


F18093-5744N 


3 


N 


Y 


N 


Y 


Y 


N 


Medium 


F18093-5744 C 


2 


N 


Y 


Y 


Y 


Y 


N 


Medium 


23128-5919 


1 


Y 


Y 


N 


7 


N 


N 


Low 


16007+3743 


Rl 


Y 


Y 


Y 


Y 


Y 


N 


High 




R2 


Y 


Y 


Y 


Y 


N 


N 


High-medium 




R3 


Y 


Y 


Y 


Y 


Y 


N 


High 



Notes. The different criteria are explained throughout the text. The letter Y indicates that the given complex satisfies the criterion, 
while the letter N indicates that it fails. Symbols with a question mark are either doubtful or indicate that we do not have the data 
to study the criterion. The last column indicates the probability that a given complex constitutes a TDG candidate based on these 
criteria. 



